Noise activated granular dynamics 



Fabio Cecconi, 1 Andrea Puglisi, 1 Umberto Marini Bettolo Marconi, 2 and Angelo Vulpiani 1 

'Dipartimento di Fisica, Universita "La Sapienza" , P.le A. Moro 2, 1-00185 Rome, Italy and 
INFM Center for Statistical Mechanics and Complexity, Italy 
2 Dipartimento di Matematica e Fisica, Universita di Camerino, 
Via Madonna delle Carceri, 1-62032 , Camerino, Italy and INFM, Unita di Camerino 

We study the behavior of two particles moving in a bistable potential, colliding inelastically 
with each other and driven by a stochastic heat bath. The system has the tendency to clusterize, 
placing the particles in the same well at low drivings, and to fill all of the available space at high 
temperatures. We show that the hopping over the potential barrier occurs following the Arrhenius 
rate, where the heat bath temperature is replaced by the granular temperature. Moreover, within 
the clusterized "phase" one encounters two different scenarios: for moderate inelasticity, the jumps 
from one well to the other involve one particle at a time, whereas for strong inelasticity the two 
particles hop simultaneously. 
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Granular gases Q, i.e. assemblies of inelastic particles 
losing a little kinetic energy at each collision, exhibit a 
variety of complex behaviors, such as clustering Q, spon- 
taneous formation of vortices Q , lack of energy equipar- 
tition 4], non-Maxwellian velocity distributions and 
so on, which provide new challenges to statistical me- 
chanics. Analogies between standard condensed matter 
and granular matter remain a valuable route to improve 
our knowledge of the latter 0. As an instance, we men- 
tion the analogy between the granular temperature, Tq, 
defined as the average kinetic energy per grain, and the 
temperature of molecular gases. It has been observed 
that the equality of the granular temperatures of two 
granular gases is not the condition for thermal equilib- 
rium between them This is not too surprising since 
it reflects the non-equilibrium nature of granular sys- 
tems. Nevertheless, one can ask whether Tq possesses 
other useful properties of the thermodynamic tempera- 
ture, such as that of controlling the rates of activated 
processes and the direction of energy fluxes. The answer 
is relevant in the construction of the hydrodynamics of 
granular systems. 

The present work is inspired by the experiment probing 
the behavior of vibrated sand in a vertical box made up 
of two identical compartments communicating through 
a small orifice located at a certain height. One observes 
that for vigorous shaking the two halves are equally popu- 
lated, whereas below a critical driving intensity the sym- 
metry is broken Q. The existing theoretical explana- 
tions are based on hydrodynamic descriptions. Accord- 
ing to Refs. @>0,^3> the behavior can be captured by a 
phenomenological mesoscopic flux model, depending on 
a parameter which is a function of the inelasticity, the 
driving intensity and the number of particles. One as- 
sumes the existence of a stationary state and, focusing 
on slow variables, neglects the role of temporal and spa- 
tial fluctuations. 

Our treatment, instead, represents a shift from the hy- 



drodynamic to the statistical mechanical level, in which 
quantities such as temperature and density fluctuations 
are obtained in terms of the microscopic coordinates of 
the particles and their interactions. In a simplified model, 
we relate the crossover from the equally populated phase 
to a broken symmetry phase to the existence of two dif- 
ferent kinds of thermally activated processes. This is a 
strong indication that the kinetic temperature is a pa- 
rameter characterizing, not only velocity fluctuations, 
but also the activation dynamics across energy barriers. 

Futhermore, the simplicity of the model allows us to 
make explicit the dependence of the kinetic temperature 
on the relevant parameters of the model. 

Our model is closely connected to that, proposed over 
sixty years ago by H.A. Kramers 0,0], describing a 
reaction occuring via a thermally activated barrier cross- 
ing. We consider two inelastic hard rods (the simplest 
granular gas) bound to move on a line in the presence 
of a bistable external potential U(x) = —ax 2 /2 + bx 4 /4:, 
(mimicking the compartmentalized box). The particles 
are coupled to a bath which exerts upon them a velocity 
dependent friction and a random force. In the absence 
of collisions the particles evolve according to: 
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where, prime indicates the spatial derivative, Xi (i = 1,2) 
represents the position of particles and 7 is a friction co- 
efficient. The Gaussian noise £i(t), with < £j(t) >= 
and < &(t)£,-(s) >= 2MK b "fT b 8ij8{t - s), describes the 
exchange of energy with the surroundings |l3j schema- 
tized by a heat-bath of temperature TJ,. In the following 
we set M = 1, K B = 1, a = 1.5, b = 0.05 and 7 = 0.1. 
The relevant parameters are: the positions of minima 
Xmin = ±-\/a/6 = L/2, the location of the maximum, 
Xmax = 0, and the energy barrier, AU = a 2 /4b = 11.25, 
together with the curvatures = U"(x m in) — 2a and 
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Upon taking into account collisions, the velocities of 
particles change instantaneously when their separation 
x 2 — x\ equals the hard core diameter, d = 1, according 
to the rule: 
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where < r < 1 is a restitution coefficient and i =/= j = 
1,2. Let us recall that in the non- interacting case (no 
collisions), each particle on the average sojourns a time 
given by 



t = Aexp(AU/T b ), 



(3) 



where the prefactor A has been calculated by many au- 
thors since the work of Kramers 0] . 

The basic phenomenology of the model is illustrated 
in Fig. n The relative distance, y = x 2 — X\, between 
particles fluctuates in time showing time intervals of av- 
erage lifetime t 2 , when they are confined to the same 
well (y ~ d) alternated with intervals, of average lifetime 
ti, when they sojourn in separate wells (y ~ L). Thus, 
the system behaviour shows the existence of two differ- 
ent time scales n and t 2 employed to characterize two 
regimes. 




FIG. 1: Relative distance xi — x\ as a function of time for a 
system with r — 0.9 and Tb = 4.0. The solid line indicates 
the diameter of the rods, the dashed marks the well separation 
L ~ 10.95. 



We shall see that, as the driving intensity, Tb, or in- 
elasticity, r, are varied, the system undergoes a crossover 
from the regime in which particles are far apart most of 
the time (i.e. t± > t 2 ), to a clusterized regime charac- 
terized by ti < T2. Moreover, we shall show that the 
dependence of T2 and T\ on the model parameters can be 
captured by a simple extension of formula (|3J), replacing 
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where t\ and t 2 are the times the particles spend, dur- 
ing [0,t], in different wells or in the same one respectively. 
The temperature T2 represents the velocity variance con- 
ditioned to the fact that the two particles belong to the 
same well, whereas T\ is the same quantity when these 
move in different wells. 

Two physical effects are present: the hard core repul- 
sion and the inelasticity of collisions. For the sake of 
clarity, let us begin the dicussion with the elastic sys- 
tem (r = 1). In this case, we deal with an equilibrium 
system - the measured T2 and T\ coincide with the heat 
bath temperature T b |14| - and therefore, we expect that 
the most probable configuration minimizes the free en- 
ergy. This configuration is constituted by a single particle 
in each well, and corresponds to the condition t\ > t 2 . 
The results of simulations, shown in Fig. verify this 
scenario. How do we quantify these escape times? As 
displayed in Fig. [3 t 2 and t\ still follow the Arrhenius 
exponential behavior of Eq. Q, however, with a suitable 
paremeter renormalization: 
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where k — (1,2) indicates single or double occupation, 
W x = AU and W 2 = AU - 6U < AU. The correc- 
tion SU takes into account the effect of the excluded 
volume repulsion: when two grains belong to the same 
well their center of mass lies higher than if they were 
in separate wells. This determines a reduction SU — 
a(d/2) 2 + b/A(d/2) 4 of the effective energy barrier and 
makes t 2 < t\. This is a typical correlation effect, be- 
cause the repulsion renders less likely, with respect to the 
non interacting case, the double occupancy of a well. The 
smaller the ratio between the well width and the particle 
diameter, the stronger the reduction of the escape time 

The above scenario changes in the presence of inelas- 
ticity (r < 1) because dissipation tends to promote the 
double occupation of a well. Thus, upon lowering the 
temperature T b , we expect a crossover from the regime, 
where the double occupancy of a well is unfavoured (i.e. 
t 2 < ti), to the regime, where particles spend most of 
the time together in the same well i.e. t 2 > t\. In Fig. [21 
for inelasticity r = 0.9, this crossover is observed and 
occurs at T b ~ 10.0 (upper inset). This behavior is the 
analogue of that reported by several authors 0, E| . The 
origin of the crossover lies on the fact that, in the inelas- 
tic system, temperatures T2 and T\ are no more equal to 
T b and furthermore T2 <T\. Thus, the mean lifetime of 
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FIG. 2: Arrhenius plot of mean escape times r. Open sym- 
bols refer to elastic case: the escape time is ri (circles) when 
a well is singly occupied, T2 (squares) when a well is doubly 
occupied. Full symbols, instead, correspond to the inelastic 
system (r = 0.9). Linear behaviour indicates the validity of 
Kramers theory with renormalized parameters and the slopes 
agree with values obtained from Eq. JSJ. Upper inset: en- 
largement of the crossover region. Lower inset: Data collapse 
of distribution of escape times for different temperatures. The 
arguments of the exponentials (dashed lines) in the same fig- 
ure have been obtained by applying formula ©. The overall 
agreement between the above prediction and the values of ti 
and T2 obtained by simulation is rather good. 



the clusterized and non clusterized regimes can be still 
described by expressions ijBJ), but now the temperature 
difference competes with the excluded volume correction, 
eventually leading to T2 > t\ . A simple argument can be 
used to estimate the shift of T2 from T&. For moderate 
driving intensity, T1 is nearly equal to Tf,, while T2 is 
lower than Tf, by a factor which depends on the inelastic- 
ity. In Fig. [3] we show these temperatures as functions of 
T},. It can be observed in Fig. [5] that T2 varies linearly 
with Tf, and its slope is a decreasing function of the in- 
elasticity (1 — r). A good estimate of temperature T2 can 
be obtained by considering the two particles in a single 
harmonic well V(x) = u>^ in x /2. Balancing the power 
dissipated by collisions —(1 — r 2 )v(v2 — ui) 2 /4 and vis- 
cous damping (— 27T2) with the power supplied by the 
external source (27T,) we arrive to the expression: 
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where the collision frequency v is estimated as uj m i n /ii 
(the factor 2 stems from the excluded volume effect). 
We have assumed that the precollisional relative velocity 
(v2 - vi) 2 ~ (< v\ > + < v\ >) = 2T 2 . An improved 
value of T2 results when using the value of v obtained 
from the simulation. The dependence of the formula on 
the restitution coefficient is shown in the inset of Fig. [21 



In the lower inset of Fig. [5] we plot the probability 
distributions of escape times t\ and T2 for several simu- 
lations. All the distributions are characterized by a peak 
at the origin and an exponential tail. When rescaled to 
have the same average, all the tails collapse to a single 
curve. Such exponential tails are typical of the original 
Kramers model for thermally activated barrier crossing. 
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FIG. 3: Plot of the temperature T2 vs. T& for various choices 
of the elasticity parameter r. From top to bottom the curves 
indicate r = 0.98, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4. As r grows toward 
the value 1 the lines approach the diagonal (solid line) where 
T2 = Tf,. The inset shows the ratio T2/T, obtained by the 
slopes of T2 vs. T, for the same values of r. The continuous 
line indicates the theoretical estimate given by formula £J. 

A new non trivial phenomenon occurs at low tempera- 
ture and small restitution coefficient. The escapes from a 
doubly occupied well become correlated, i.e. the transi- 
tion takes place as a collective motion of the two particles. 
In other words, when the inelasticity is strong the relative 
motion of the two particles, due to repeated collisions, 
becomes frozen and they tend to form a "molecule". It 
is easy to show (within an harmonic treatment of the 
potential) that the noise acting on the center of mass co- 
ordinate corresponds to a reduced heat bath temperature 
Tf,/2. We computed the exit time r m for the "molecule" 
and compared it with the characteristic times r 2 for sev- 
eral values of the inelasticity. As shown in Fig. 0] when 
the driving temperature decreases, t 2 corresponding to a 
strongly inelastic system (r = 0.5) exceeds r m . This is 
the signature that the most probable evolution of a dou- 
bly occupied state involves the simultaneous hopping of 
the two particles in the adjacent well, without breaking 
the pair. This freezing of the internal motion, is related 
to the problem of the dynamics of two randomly accel- 
erated particles on a line considerd in Refs. [lfil Ht| |. or 
equivalently, to the problem of a single particle moving 
on the half line (x > 0) in the presence of an inelastic 
(impenetrable) wall. The authors predicted that at fixed 
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driving intensity and for r below a critical value, r = r c , 
the particle localises at the wall. We argue that the lo- 
calization mechanism is identical to that leading to the 
formation of the "molecule" in our model. 




FIG. 4: Arrhenius plot of the mean escape times T2 versus 
1/T& for several values of r. Thick solid line indicates the 
escape time r m of single particle of mass 2M corresponding to 
simultaneous jumps. At about Tb = 10 and r — 0.5 we observe 
a crossover between T2 and r m signalling the formation of a 
stable molecule. The inset shows the difference r m — T2 vs. 
Tb to pinpoint the crossover where 

In summary, whereas previous studies on comparti- 
mentalized systems were based on coarse grained de- 
scriptions of granular gases, our approach focuses on the 
statistical description of particle motions by means of a 
characterization of their stochastic fluctuations. The ba- 
sic phenomenology illustrated above, suggests that after 
a suitable redefinition of parameters (T b replaced by T 2 
and AU by W%) the present problem can be mapped onto 
that of a single particle hopping between two wells at dif- 
ferent temperatures T\ and T 2 . We have shown how T 2 
decreases with the inelasticity. The local granular tem- 
perature is the relevant control parameter in determining 
the direction of the energy flux (particles flow from hot 
places to cold places). The activation rates, t 2 _1 and r-f 1 
display an Arrhenius dependence on the temperatures Ti 
and T 2 , respectively. However, as r decreases further, a 
deeper scrutiny reveals a more complex scenario. The 
model displays a further "transition" at even low tem- 
peratures because the two rods form a "bound" pair as a 



result of extremely frequent collisions. The escape from 
the wells can occur only "in tandem" . 

As a natural extension of the present model, due to 
the continuous interest on transport properties on rough 
surfaces, compartimentalized systems and ratchets, one 
could examine the dynamics of an ensemble of grains in 
a larger number of wells. This would allow to determine 
explicitly the variation of the granular temperature and 
of the transition rates with the occupation number of the 
wells. The latter are key ingredients in hydrodynamic 
descriptions. 
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